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We consider the Jaynes-Cummings model describing the interaction of a Cooper pair box (GPB) 
and a nanoresonator (NR) in the presence of a Kerr medium and losses The evolution of the entropy 
of both subsystems and the GPB population inversion were calculated numerically. It is found that 
these properties increase when the NR frequency is time-dependent, even in the presence of losses; 
the effect is very sensitive to detuning and disappears in the resonant regime. The roles played by 
the losses affecting the GPB and the NR are also compared. 
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I. INTRODUCTION 

The implementation of the quantum computer is con¬ 
sidered as a big challenge in the recent years due to the 
failure of classical computers ability to simulate and com¬ 
pute the large scale data. The interaction of the sin¬ 
gle photon field with superconducting circuit to generate 
quantum bits is studied experimentally mia and the¬ 
oretically mg. Wallraff et al. have observed a strong 
coupling between the superconducting two-level system 
and a single microwave photon which prove the gener¬ 
ation of qubits. Recently, the construction of the first 
quantum computer by Canadian company (D-Wave) us¬ 
ing the system of the superconducting qubits was shown. 
This quantum computer succeeds to solve famous prob¬ 
lem in mathematics, the Ramsey numbers [6]. 

The study of the interaction between a two-level atom 
and a radiation field is very important in quantum me¬ 
chanics, as its applications in laser physics and quantum 
optics. This system has similarities with others, e.g., 
the interaction between a Cooper pair box {CPB) and a 
nanomechanical resonator {NR). There are many works 
in the literature dealing with these systems [71 IHl, but 
only few of them treat these systems in the situation 
where one of the frequencies U HQ], or the amplitude 
m, varies with time. In the CPB-NR system these 
variations change the subsystems coupling and modify 
their dynamic properties, e.g., showing amplification of 
the excitation transition rate in subsystems m- 

An important focus of quantum optics is concerned 
with the atom-field coupled system. Inspired by vari¬ 
ous tests applied to this system, and its limitations, the 
researchers have passed from the light domain to the mi¬ 
crowave domain of superconducting cavities coupled to 
Rydberg atoms, or to the quantum electrodynamics cir- 
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cuit coupled to nanoresonators. After being extended 
to broader context the system was used to investigate 
Landau-Zener transitions m-, atomic physics and quan¬ 
tum optics [14]; mechanisms for photon generation from 
quantum vacuum m-, quantum simulation, challenges, 
and promises of fast-growing field m- Here we will as¬ 
sume the laboratory substituting the atom (field) by the 
CPB {NR) [17]. This system have been analyzed in sev¬ 
eral works, e.g.: quantum network m-, phonon blockade 
m-, squeezed states and entangled states [20]; cooling 
mechanical oscillators [21] ; Bell inequality violations [22] . 

In this report we will employ the Jaynes-Cummings 
model to treat a CPB-NR system with losses; a nonlin¬ 
ear Kerr medium is also added to include influences of 
nonlinearities. We investigate the time evolution of the 
CPB population inversion {I{t)), as well as the statisti¬ 
cal properties of both subsystems. We consider the NR 
initially in a coherent state and the CPB in its excited 
state. The influence caused by the time-dependent NR 
frequency upon the properties of both subsystems is stud¬ 
ied. We consider the combined effects of nonlinearity and 
losses upon the dynamics of the CPB population inver¬ 
sion and upon the NR entropy, the latter being used as 
a measure for the degree of entanglement. We compare 
the results obtained in the resonant case {ujnr = ^cpb) 
with those obtained for small detunings between the CPB 
and NR. The influences of losses from the NR and CPB 
upon the mentioned properties are also compared. 

II. THE HAMILTONIAN SYSTEM 

A superconductor CPB charge qubit is adjusted to the 
input voltage Vg of the system through a capacitor with 
an input capacitance Ci. Observing the configuration 
shown in Fig. 0 we see three loops: a small loop in 
the left, another in the right, and a great loop in the 
center. The control of the external parameters of the 
system can be implemented via the input voltage Vg and 
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the three external fluxes and $e- The control of 

theses parameters allows us to make the coupling between 
the CPB and the NR. In this way one can induce small 
neighboring loops. The great loop contains the NR and 
its effective area in the center of the apparatus changes 
as the NR oscillates, which creates an external flux ^e(^) 
that provides the coupling between the CPB and the NR. 


We will consider H = 1 and assume as identical the 
four Josephson junctions, with the same Josephson en¬ 
ergy Ej, the same being assumed for the external fluxes 
and i.e., with same magnitude, but opposite sign: 

la this way we can write the Hamilto¬ 
nian describing the entire system as 


H = ujQa^a -h 4:Ec 



&z — 4^j cos 



cos 



( 1 ) 


where d) (a) is the creation (annihilation) operator for 
the NR excitations, corresponding to the frequency uj 
and mass m; E^j and E^ are respectively the energy of 
each Josephson junction and the charge energy of a single 
electron; Ci and Cj stand for the input capacitance and 
the capacitance of each Josephson tunnel, respectively; 
^0 = ^ is the quantum flux and Ng = is the charge 
number in the input with the input voltage Vg. We have 
used the Pauli matrices to describe our system operators, 
where the states | 0 ) and | 1 ) represent the number of extra 
Cooper pairs in the superconducting island. So we have, 
= | 1 )( 1 | - | 0 )( 0 | and • 

The magnetic flux can be written as the sum of two 


terms, 

= 4>6 + Blx , ( 2 ) 

where the term ^5 is the induced flux, corresponding to 
the equilibrium position of the NR and the second term 
describes the contribution due to the NR vibration; B 
represents the magnetic field created in the loop and £ 
is the length of the NR. We write the displacement x 
as X = xo{a'^ + a), where xq represents the NR ampli¬ 
tude oscillation. Substituting the Eq. (§ in Eq. § and 

controlling the flux ^5 we can adjust cos = 0 to 

obtain 


H = ujQCL^a + 4Ec 



and making the approximation << 1 we find the 

Hamiltonian in the form, H = ujod'^a -h ^ujc^z + Ao(a't' + 
d)dx^ where the constant coupling Aq and the effective 

energy Uc are given by Aq = cos (^) , 

ujc = 8 ^c “ I) • rotating wave approximation 

the above Hamiltonian results in the form, 

Heff = + ^ujcdz + Ao((J+d + (4) 


— 4E j cos 



/ 'kB£x\ 


(3) 


Next, we consider a more general scenario, substituting 
CU = LUq E f (t) , Xq ^ A(t) = Ao(l + 

and x(t) = XoP^f{t) [23l|24]. In addition we consider the 
presence of the term K{t) standing for the time-dependent 
loss affecting the CPB, the term 6{t) being the same 
for the NR, and x(t) is the response time of the Kerr 
medium. This extended and somewhat realistic scenario 
requires the substitution of the Hamiltonian Hef f by the 
Hamiltonian Jif given by 


= uj{t)d^d -h ^ujc{t)d-z + A(t) (d(j+ + -h “ ^S{t)d^d. (5) 


In the Eq. (§ the first (second) term describes the NR 
{CPB) subsystem, the third term represents the CPB- 
NR coupling, the fourth term stands for the time de¬ 
pendence of the Kerr medium, the fifth and the sixth 
terms, n{t) and 5{t), were mentioned above. The oper¬ 
ators (T± = (Tx E icFy leads to (j+ = |1)(0|, d_ = |0)(1|; 


n{t) is the CPB decay coefficient from the excited level 
| 1 ) to the fundamental state | 0 ), 5{t) is the same for the 
NR and X{t) is the time-dependent coupling between the 
CPB and NR. Actually, the inclusion of loss in the system 
turns its treatment somewhat more realistic since dissi¬ 
pation is ubiquitous in the real world; as consequence 
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FIG. 1. Model for the CPB-NR eoupling. 


the Eq. corresponds to the evolution of a system 
described by a non-Hermitian Hamiltonian [25H29] . On 
the other hand, as the response time of the Kerr medium 
is assumed so fast the medium follows the NR adiabati- 
cally and the third-order nonlinear susceptibility can be 
modulated by the NR frequency uj {t). 

The wave function that describes our system can be 
written as, 


oo 

1^ (t)) = E \^,n) + Co,n (t) | 0 , n)], ( 6 ) 

n 

where Ci^n {t) and (7o,n {t) are respectively the probabil¬ 
ity amplitudes of the states |l,n) and | 0 ,n), namely, the 
CPB in its excited state |1) or ground state |0) with el¬ 
ementary excitations n in the NR. As mentioned before, 
at t = 0 the system is decoupled, the CPB initially in 
its excited state |1) and the NR in a coherent state ja). 


given by 


i«) = E • 

n=0 


(7) 


The total wave function in the initial state can be 
written as |T( 0 )) = ? with the initial am- 

plitudes C'o,„(0) = 0 and = 1 - The 

Schrodinger equation for the present system, described 
by non-Hermitian and time-dependent Hamiltonian in 
the Eq. <§ is. 


dt 




(8) 


We obtain the following set of coupled equations 
of motion for the amplitude probabilities Ci^n (t) and 
Co^n-\-l (^)* 


^ - n) - + nS{t)fj Ci^n {t) - iX{t)^{n + l)Co,n+i (t), (9) 


dCo,n+l (t) 


= ( -i{n + + 




■ n) - (n- 


l)d(i) ) Co,„+i (t) - iX{t)y'{n+ l)Ci,„ (t), (10) 


whose solutions allow us to calculate the entropy of the III. TIME EVOLUTION OF THE POPULATION 
NR subsystem and the population inversion of the CPB. INVERSION 


The present approach also allows us to investigate the 
CPB dynamics in a non perturbative way. A convenient 
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way to characterize the response to the NR influence is 
given by the CPB population inversion. This parameter 
is defined as 


n=0 


( 11 ) 


To calculate this property I{t) we will assume the 
NR frequency varying with time as uj {t) = c^o + f{t). 
The third order nonlinear susceptibility is modulated as 
x{t) = Xo + ^f{t) and we also assume the NR initially in 
a coherent state with the mean excitation number n = 25 
and ^ ^ = 20k. We consider the time evolution of 

the population inversion for different values of the decay 
coeficients K{t) and S{t). 


(a) 



A. Resonant case: /(t) = 0 

If the model is restricted to the usual case, with the NR 
frequency constant, the Eqs. ([^and[TQ|) can be solved nu¬ 
merically. In Figure ([^a) we plotted the inversion I{t) as 
a function of time. In absence of losses the average inver¬ 
sion, defined by its value during the collapse, is greater 
than zero, close to 0.5, when we consider the value ^ 
= 0.2. However, as well known, neglecting the presence 
of the Kerr medium = 0), the value of the average 
inversion vanishes. 

Here the third-order nonlinear susceptibility represents 
the coupling of the NR and the Kerr medium. The higher 
the value of x the stronger the coupling of the NR with 
the Kerr medium, and reversely. Considering the case 
^ = 0.2 and the system with no loss (/^(t) = S{t) = 0 ) we 
observe the collapse-revival effect (see Fig.([^a). In this 
figure, the horizontal line crossing the value I{t) = 0.5 of 
the average inversion is due to the presence of nonlinear¬ 
ity, whose absence turns the average inversion null. We 
note in this figure the amplitude of oscillation decreas¬ 
ing with time, due to the presence of the Kerr medium. 
However, when considering only the CPB loss (^c(t) ^ 0, 
5{t) = 0 ) the collapse-revival effect reappears, with sub¬ 
sequent oscillations that smoothly vanish (see Fig. ^ 
b)). Contrarily, when considering only loss in the NR 
{K,{t) = 0 , S{t) 7 ^ 0 ) these oscillations rapidly vanish. So 
the spoiling effect caused by the NR loss dominates that 
coming from the CPB loss, see Fig. ([^c) and (j^b). 
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B. OfF-resonant case: /(t) = r sin(ct;/) 

When the system is non resonant, with ^ = 10 and 
^ = 1, and assuming only CPB loss 7 ^ 0, 5{t) = 0), 
we observe periodical behavior of the population inver¬ 
sion, its amplitude decaying over time - see Fig.(|^a); in 
this case the collapse-revival effect is not observed. How¬ 
ever, this behavior is modified when the loss comes from 
the NR\ in this case the inversion and the oscillations 
are destroyed. Considering uj' increasing and only the 


FIG. 2. Time evolution of the population inversion for dif¬ 
ferent values of the parameters Hi{t) and 5{t) for: (n) = 25, 
ct;o/Ao = ooc/N = 20k, f{t) = 0; Xo/Ao = 0.2; (a) k/Xq = 0.0 
and S/Xq = 0.0; (b) k/Xq = 0.01 and (5/Ao = 0.0; (c) 
k/Xq = 0.0 and S/Xq = 0.01. 


CPB loss - see Fig. ib), the system continues displaying 
the population inversion and oscillations, the latter ex¬ 
hibiting small periods, an effect not shown when the loss 
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comes only from the NR - see Fig.(|^c). 


IV. TIME EVOLUTION OF THE ENTROPY 



Nowdays some devices are based on quantum mechani¬ 
cal phenomena, and this holds also for information trans¬ 
mission. For example, in optical communication a polar¬ 
ized photon can carry information. Now, the entropy of 
entanglement is defined for pure states as the von Neu¬ 
mann entropy of one of the reduced states, e.g., the NR 
entropy, 

SNR{t) = - , 

( 12 ) 

where = 1{{C\C) + {S\S) ± 5 [((C|C) - 

(5|5))2 +4|(C|5)|^]V2) with {C\C) = E:Lo \CiAt)\\ 
{S\S) = and (C|5) = (5|C)* = 

TA=^Cln+l{t)A,n+l{t). 

In our calculations using the Eq. we will assume 
the NR frequency varying in the form uj (t) = c^o + /(t). 
The third order nonlinear susceptibility is modulated in 
the form x(t) = Xo + ^/(^) we also assume the ini¬ 
tial NR in a coherent state with the mean number of 
photon n = 25 and ^ ^ = 20/c. Next we will con¬ 

sider the time evolution of the entropy for different values 
of the decay coefficients and d(t). 


A. Resonant case: f{t) — 0 

We will use the values of parameters in Fig. In 

an ideal system the parameters n{t) and 5{t) are null, as 
assumed in Fig. (ga); in this figure the maximum of the 
NR entropy is close to In 2; after the start of the inter¬ 
action the NR entropy gradually tends to its minimum, 
then returns to its maximum and remains oscillating reg¬ 
ularly. 

For small values of decay in the CPB^ as ^ = 0.01, 
and ideal NR {S = 0) for the case ^ = 0 . 2 , the maxi¬ 
mum value of the entropy shows no significative changes 
for small times - see Fig. (|^ b). For larger values of 
time the amplitude of the entropy oscillations decreases, 
maintaining the periodicity. If instead we add a small loss 
only in the NR^ say ^ = 0 . 01 , the entropy oscillations 
vanishes rapidly (cf. Fig. (|^c)), showing the entropy 
being more sensitive to the loss in the NR than that in 
the CPB. For larger value of 6{t) and hz{t) the entropy 
movies rapidly to zero, as expected, due to the passage 
of both subsystems to their respective ground states. 


B. OfF-resonant case: f(t) =Tsin(ct;/) 


FIG. 3. Time evolution of the population inversion for dif¬ 
ferent values of the parameters K,{t) and S{t) for: (n) = 25, 
^o/N = ujc/Xo = 20A:, s/Xq = 0.001; t/Xq = 10; f{t) = 
Tsin(cc;/t); Xo/Ao = 0.2; (a) k/Xq = 0.01; S/Xq = 0.0 and 
ujf/Xo = 1; (b) k/Xo = 0.01; (5/Ao = 0.0 and ujf/Xo = 20; (c) 
k/Xq = 0.0; 6/Xo = 0.01 and uf/Xo = 20. 


Let us now consider the variation in the detuning pa¬ 
rameter, where r and u;/ are parameters that modulates 
the NR frequency. Our discussion is limited to the con¬ 
dition r ujc^ ^0 and also assuming that uof is small 
to avoid interaction of the CPB with other modes of the 
NR. We have chosen various values r of amplitude modu¬ 
lations to verify the entanglement properties between the 



















































































6 



FIG. 4. Time evolution of the entropy for different values of 
the parameters and S{t) for: (n) = 25, Uc/^o = coq/Xq = 
20/c, Xo/Ao = 0.2, /(t) = 0. (a) k/Xq = 0,0Ao and S/Xq = 
0.0; (b) k/Xq = 0.01 and (5/Ao = 0.0; (c) k/Xq = 0.0 and 
S/Xo = 0.01. 


FIG. 5. Time evolution of the entropy for different values of 
the parameters K(t) and S{t) for: (n) = 25, lUc/Xq = loq/Xq = 
20/c, Xo/Ao = 0.2, r/Ao = 10, and e/Xo = 0.001. (a) k/Xq = 
0.01, 6/Xo = 0.0, and ujf/Xo = 1; (b) k/Xq = 0.01, (5/Ao = 
0.0, and uf/Xo = 20; (c) k/Xq = 0.0, (5/Ao = 0.01, and uf/Xo = 
20. 


CPB and NR. We also use various values cj/ of frequency 
modulation to see its influence upon the CPB-NR entan¬ 
glement - see Figs. (§. Comparing the Fig. (|^a) with 


the Fig. ([^a), we note that the entropy loses its peri¬ 
odical oscillations, and when u;/ increases this behavior 
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becomes more evident - see Fig. §b). Now, comparing 
the entropy in Fig. (|^c) with that in Fig. ([^c), keeping 
the same parameters, we see the entropy going to zero 
in both cases, but the Fig. ^c) displays an interesting 
effect: even in presence losses, the maximum value of the 
entropy grows, reaching the value In 2, after which it os¬ 
cillates downward. This effect came from the sinusoidal 
modulation of the NR frequency 

V. CONCLUSION 

In the present work we have considered the interaction 
of a CPB and an NR in the presence of a Kerr medium 
and losses affecting both subsystems. Concerning the 
inffuence of the loss affecting the entropy of both subsys¬ 
tems and the population inversion of the CPB we have 
observed the dominant role played by the NR upon that 
played by the CPB. The dissipation causes deterioration 
of the CPB excited level whereas convenient modulations 
favors the control of certain properties of the system. It 


was also shown that certain choice of the time-dependent 
frequency makes higher the maximum value of the en¬ 
tropy, even in the presence of dissipation (see Figs. (|^c) 
and Fig. ^c) ) the same occurring for the CPB popula¬ 
tion inversion - see Figs, ([^c) and Fig. (j^c). Concerning 
the entropy, this result is very important for information 
transmission since the transmission of maximum infor¬ 
mation through a quantum channel is exactly the von 
Neumann entropy [301EB. These effects are very sen¬ 
sitive to detuning 0.04%) and disappear in resonant 
regime. The results suggest that it is possible to perform 
a dynamic control of certain properties of this system, 
via convenient manipulation of the parameters. We hope 
that these results might shed light in this scenario, fur¬ 
nishing new insights for researchers in this area. 
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